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By Max A. Heaslet and Otway O’M. Pardee 


SUMMARY 


The critical Mach number, as a function of lift coef- 
ficient, is determined for certain thin and moderately thick 
NACA low— drag airfoils. The results, which are given graphi- 
cally, include calculations on the same airfoil sections with 
plain flaps for small flap deflections. Curves are presented 
indicating optimum critical conditions for the airfoils with 
flaps and are in a form so that they may he compared with 
corresponding results for zero flap deflections. 

The calculations indicate that, through the use of plain 
flaps, an increase may he realized in the 1 if t— co ef f ic i ent 
range for which the critical Mach number is in the region of 
high values characteristic of low— drag airfoils. 


INTRODUCTION 


The necessity of attaining higher speeds and higher 
altitudes in military aircraft has fecused increasing atten- 
tion on the critical speeds of the airfoil sections used and, 
as a result of investigations concerned with the calculation 
of these critical speeds, certain properties of favorable 
airfoil sections hare become known. For example, it is 
possible to say as a general conclusion that the type of 
pressure distribution associated with low— drag airfoils is 
one which is also favorable to the production of h igh— cr it i cal 
speed characteristics. Moreover, as pointed out in refer- 
ence 1, the thickness ratio and the camber of a wing section 
play an important role since the maximum value of critical 
Mach number decreases approximately linearly as the camber 
and thickness of an airfoil increases. It has therefore been 
almost inevitable that the airfoil sections used on recently 
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designed high-speed aircraft have had the type of pressure 
distribution associated with the low— drag airfoil and have 
been thinner than those used formerly. 

Or it ical— speed curves, such as are given in reference 1, 
show that the low-drag airfoil has a region of lift coef- 
ficients, more or less symmetrically disposed with respect 
to the design lift coefficient, in which the critical Mach 
number variation is small, the maximum critical Mach number 
being achieved within the region. Outside this sector, 
which corresponds roughly to the 1 if t— coef f ic i ent range for 
which the low— drag properties of the airfoil hold, there is 
a sharp decrease in the value of critical Mach number. It is 
obvious that a particularly advantageous situation exists if 
it is possible to design the wing section of an airplane so 
that the h igh—cr it ical— speed and low-drag regions of the 
wing extend beyond the 1 if t— coef f icient range for normal 
operations. The difficulty of achieving this has already 
been encountered in the design of fighter aircraft where 
demands on the maneuverability at high speeds are great. 

The problem arises again in the case of long-range bombers 
since, on extended flights with attendant fuel consumption 
and with the accompanying disposal of bomb loads, the vari- 
ation of lift coefficients required may be quite large. The 
situation is particularly acute for Jet-propelled bombers 
since high speeds are possible of attainment over a wide range 
of altitude. 

The theoretical results of reference 1 show quite clearly 
that, as the thickness of an airfoil section decreases, the 
maximum critical Mach number is increased but that this in- 
crease is brought about at the expense of the 1 if t— co ef f i c lent 
range for the h igh—cr it ical— speed region. Thus, the high- 
speed requirements in the design of an airplane may call for 
a thin wing section while other specifications may be such 
that the extent of lift coefficients needed at high speeds 
extends beyond the natural range of the airfoil. As a conse- 
quence, it becomes highly desirable to investigate any method 
whereby an extension of this range may be effected. One such 
method which could presumably be used for this purpose is 
the use of full— span plain flaps, for in this manner the 
camber may be modified and the load distribution over the air- 
foil disposed so that the sharp growth of the pressure peak 
near the nose is restricted. 

In the present report, calculations have been carried 
out at the request of the Air Technical Service Command, 

U. S. Army Air Forces, to determine the critical Mach numbers. 
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as a function of lift coefficients for the NACA 64—, 65—, 
and 66— series low— drag airfoil sections with th ickness— chor d 
ratios equal to 0.06, 0.08, 0.10, and 0.13. These sections 
have constant ideal lift coefficients of 0.3 at which lift 
the load is distributed uniformly over the chords. Another 
portion of the theoretical calculations 18 devoted to the 
determination of the critical Mach numbers of these same 
airfoil sections with plain flaps, for small flap deflections, 
in an attempt to study the effect of such flaps on the criti- 
cal Mach number curves of the airfoils. 

An experimental investigation has also been carried out, 
under the same general research program, to determine the 
Mach numbers at which the force and moment characteristics 
of the same sections with flaps are divergent. This investiga- 
tion is to be reported separately and will contain a compari- 
son of the theoretical critical results with the experimentally 
evaluated divergence Mach numbers. 

A complete list of symbols, as used throughout this report, 
may be found in the appendix. 

ANALYS IS 

Computation of Critical Mach Number 


It Is now an established convention to define the critical 
Mach number of a body as the Mach number of the free stream 
for which, at some point on the surface of the body, the 
fluid first reaches a velocity equal to the local velocity 
of sound. In an analogous manner the critical compressi- 
bility speed is defined as the free— stream speed corresponding 
to that at which the critical Mach number is attained. Ex- 
perimental evidence, ebtained from the study of airfoil 
sections, indicates that compression shocks are formed 
locally on an airfoil surface soon after, if not coinci- 
dental with, the attainment of Mach numbers corresponding to 
the critical speed. This shock, however, is not well defined 
and, so far as can be observed, no strongly developed shock 
front exists until the free— stream Mach number has risen some- 
what above its critical value. It thus seems quite reason- 
able to assume, and this has been further substantiated by 
experiment, that for airfoils of limited thickness the 
critical Mach number furnishes a conservative approxi- 
mation, for the designer, for the occurrence of the flow 
breakdown which is associated with super— cr it i cal speeds and 
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produces the sudden changes in the airfoil characteristics 
that are, in general, inimical to good airplane control and 
perf orman ce. 


The critical Mach numbers in this report are calculated 
in the manner used in reference 1. Under the assumption 
that the flow is isentropic, the critical pressure coef— 
ficient i 8 given by the relation 

Y 


Per Po 


YM" 


cr 


[7 ■ 

+ k? 

y- 1 , 1 

(A Y+l 

Y+l r 

j - ij 


( 1 ) 


wher e 

p 0 static pressure in the free stream 

p cr pressure corresponding to sonic velocity at M pr 
and occurring at minimum pressure point 


V D velocity of the free stream 

a Q velocity of sound in the free stream 

M cr critical Mach number, equal to V 0 /a Q a t critical 
conditions 

Y ratio of specific heats (cp/c v = 1.4) 

p Q density of the fluid in the free stream 
dynamic pressure ( — p V 2 ^ 

0 V 2 ° 0 / 


In order that 
coefficient Pfj-o 



be related 
it is 

M«0 


to the low— speed pressure 
necessary to express the 


left-hand side of equation (l) in terms of Py iac Q an<i to 

end the Karman— Tsien formula (reference ?.) has been used. 
Equating this result, evaluated at the critical value of K , 
to the right-hand side of equation (l), gives the requisite 
expression 
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For an airfoil, or an airfoil with flap, if it is 
possible to determine the pressure distribution for a given 
low— speed lift coefficient, then the critical Mach number 
of the airfoil can be found by means of equation (2) together 
with the value of pressure coefficient at the minimum pressure 
point. As a result of such calculations, the critical Mach 
number of an airfoil is found as a function of the low— speed 
section lift coefficient, c Im = o* Por < * es *& n « however, it 

is highly desirable to know the actual lift coefficient 

corresponding to the Mach number of flight. In reference 3, 
Glauert has developed the approximation 


c l 


M 


c 4l = 0 
yi-M " 1 2 


which does relate the low— and high-speed lift coefficients, 
and in this report the G-lauert correction has been applied. 
Experimental observations show that the accuracy of this 
correction is good for Mach numbers up through critical values. 


Calculation of Pressure Distributions 
for Airfoils with Flaps 

The theoretical calculation of airfoil pressure distribu- 
tions has been outlined in reference 4, and in this reference 
tabular data are given whereby the pressure distributions may 
be calculated immediately for all standard NACA airfoil sec- 
tions with various types of camber lines. The velocity 
distribution over the airfoil is considered, in conformity 
with present theory, to be formed from three separate and 
independent parts: 

1. That part of the velocity distribution associated 

with the basic thickness form set at zero angle 
of attack 
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2. That part of the distribution associated with the 

design load of the camber line 

3. That part of the distribution associated with the 

additional load distribution and related to the 
angle of attack of the airfoil 


As a result of this 
the pressure coefficient 


theory it is possible to express 
Pj 4 = 0 in the form 


P M *0 = 1 


-( 


z ± H ± flla 


y 


(3) 


v Au Av a 

where and 

v o v o ¥ o 


are velocity ratios corresponding 


respectively to parts 1, 2, and 3. This method has been used 
throughout the present report for airfoils with and without 
flaps, and the airfoil data in reference 4 have been used in 
all cases. 


For an airfoil with plain flap it is necessary to find 

Au Av a 

the effect of the flap on the values of — and - - --- , 

o a 

and this can be best achieved by first calculating the change 
in the load distribution over the airfoil which is brought 
about by the flap deflection. In reference 5 this problem 
has been treated in a sem iempir ical fashion for the case of 
conventional airfoil sections, but the theory is not immedi- 
ately applicable to low— drag airfoils and, for this reason, 
a different approach is made modeled on the work of Glauert 
in reference 6 and Allen in reference 7. 


It is an accepted practice to divide the chordwise lift 
distribution P of an airfoil into two parts: (a) the so- 

called "basic" lift distribution P^, which depends on camber- 

line shape and is independent of the angle of attack; and (b) 
the additional lift distribution P ft , which is variable with 

angle of attack and in form is independent of the camber— line 
shape. When the flap on an airfoil is deflected, the change 
in lift distribution is called the incremental lift distri- 
bution Pg, and the two component parts are respectively, 
incremental basic distribution P^g , and incremental 
additional distribution P a g . It can be shown that the 
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incremental additional distribution due to the deflection of 
the flap is identical in form with the additional distri- 
bution for the airfoil with flaps neutral. The incremental 
basic distribution must be evaluated, however, from a knowl- 
edge of the airfoil section, the nature of the flap, and the 
flap deflection. The determination of this variation is 
therefore undertaken in the following paragraphs. 

In conformity with the assumptions usually made in 
th in—airf o il theory, the airfoil is replaced by its mean 
camber line. Figure 1 shows the assumed camber line distri- 
bution produced by the deflection of the flap, the hinge point 
of the flap lying between Xj and xg , AF approximating 

the chord— line of the airfoil. If x is measured from 
point A along AF and y is measured from A along a line normal 
to AF, then in the figure 

< < 

is linear with equation y=0 for 0 = x = Xj 

s 

is parabolic with equation y - ax +bx+d for 

< < 

Xi a X * Xg 

is linear with equation y— y a = —tan 6 ( x— xg) for 
< < 

X £ = X = c 

where a, b, and d are arbitrary coefficients and S 
is the angle the flap is deflected. The parabolic section 
is to extend over a very small portion of the camber line, 
the extent of this section being determined later, and is 
introduced to avoid the sharp break in slope which theoretically 
would exist and the subsequent requirement of a singularity 
in the velocity and lift distribution at the hinge point. 

This small portion between B and C may be thought of as a 
fairing of the camber line at the point of the sharp break 
and is consistent not only with the existence of a boundary 
layer on the surface of the airfoil, for the layer has a 
tendency to iron out such abrupt irregularities, but also with 
the geometry of the median line between the upper and lower 
surfaces of the airfoil. 


AB 

BC 

CD 
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If the transformation 

i 


X = — C ( 1 — COB 0) 
2 


is introduced, the expressions for the slopes of the three 
sections may he written in the forms 


AB: 

£ = 0 
dx 

far 0 "= 0 = 0 j 


BC: 

= 2ax + b for Xi = 

dx 

< 

x = x s 


= G 

cos0+H for 0j 

* 9 * «a 

CD: 

1 

ii 

sis 

tan 6 



= — 

< < 

A for 9 s = 8 = 

TT 

where G, H, 
are defined by 

and A , 
the above 

are introduced for 
equat ions . 

s impl i c ity 


(4) 


Reference 7 establishes the relationship 


/o p bN 

c — ) 


= _i } iz ( eot if 

2rr o dx 2 


cot 


o-e 


) 


d0 


(5) 


where 

(4 s ) 


e 

% 

dx 


basic distribution for infinitesimally thin airfoil 
at chordwise station corresponding t* 0 = 0 Q 

variable of integration 

value of 9 at an arbitrary fixed point 
slope ef camber line 


It is ebvious that this expression can be used in conjunction 
with the slepes given in equations (4) to determine the 
incremental basic distribution due to the flap deflection. 
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The integration of the integral is straightforward, and the 
f inal result is 




sin 9 o + 


A 

TT 


In 


s in-g( 9 2 + 9 0 ) 

8 ® o “ 8 J 



cos 


9 + H) 

0 



3 in-*;( 9 2^*9 0 ) s in-g( 9 9 Q ) 

sin-|(9i— 0 Q ) sini(9 s +9 0 ) 
3 3 


(6 


For a given flap deflection of 8 degrees the value of 
A is known and it merely remains to fix the coefficients 
Or and H and determine 9 t and d 2 so that the theoretical 

incremental basic distribution is consistent with experiment. 
Impose now the condition that at x = Xj (9=0 x ) the parabola 
has zero slope and a radius of curvature equal to r. This 
requires that 

c c 

q. s and H = — css 9,. 

2r 2r * 

Moreover, at x = x E ( 9 = 9 2 ) the slope of the parabola must 
equal —A , thus 

c . 

— (cos 9 2 — cos Bj) * —A 
2r 

r 

and if — , A, and 9 X are known, it is possible to find 
° 6 j + 9 2 

6a* Letting the hinge point of the flap be at — 2 — , the 

parabola can be oriented so that the values of 9 at its 
end points are symmetrically disposed with respect to 9 at 
the hinge point and for small flap deflections 9j and 8 S 
can be found. 

Using equation (6) and the derived values of the various 
parameters, a comparison was made between calculated values 
of incremental basic distribution and available experimental 
data for small flap deflections. It was found that when r 
was set equal to the thickness of the airfoil at the hinge 
point the agreement was quite good over the entire airfoil 
surface and that at the hinge point, where maximum values 
of are attained, the results were reasonably accurate 
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to Justify the use of the theory to calculate critical Mach 

numbers. 

By the methods of reference 7 it follows that 


C l b6 = 2X 8in *3 


and 

° l a6 = 2 X (tt-B b ) 


G 

+ C(0 2 — Bi)+ — (sin 
2 

+ 2H( sin 9 2 — sin 
— 2H( 9 2~0 i ) -2G( sin 


20s— sin 20!) 


e i) 

(7) 

0 2 — sin 

(8) 


where 

C[ t incremental basic lift coefficient 
l b6 

ci . incremental additional lift coefficient 
afi 

Since 0 9 i is small, for small deflections of the flap, 

it is possible to approximate equations (?) and (8) by simpler 
expressions and this was done in the calculations. Under the 
same assumption, the peak point of the Incremental basic lift 
distribution is at the hinge point and can be approximated 
quickly . 


The relation between the velocity distribution and the 
chordwise lift distribution over the airfoil has been given 
by Allen (reference 8) in the form 


and 



( 9 ) 


( 10 ) 
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where 



velocity over upper surface of airfoil 


velocity ever lower surface of airfoil 


pressure distribution over base profile 


P 


load distribution over airfoil 


Since the contribution of the flap deflection to P has been 
calculated, the effect on the velocity distribution can be 
determined. Substitution in equation (3) will give the low- 
speed pressure coefficient at any point along the airfoil; 
and, from the minimum pressure, the critical speed of the 
section is determinable. In calculating the contribution 
of the incremental basic lift distribution to the velocity 
•n the surface of the airfoil, it is well to bear in mind that 
the expression for 0 ^bS given in equation (6) was derived 

under the assumptions of thin airfoil theory. In reference 
7 it has been pointed .out .that to a first order of approximation 


p bS - <^8 


Hence, it follows that the incremental velocity associated 
with this portion of the lift distribution is given by 

—o^bS • The final form of equation (3) may therefore be 
4 

written as 


P M=0 = X 




Au ± Av a ± 


4 o F bs) 


Au 


where is determined from the design load of the basic 

V o 

camber line and — P-^g a function of the flap— chord ratio 


and the deflection angle of the flap. 
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DISCUSSION OF RESULTS 


In figure 2 are shewn the critical curves of a typical 
airfoil section (NACA 65! — 210) for various flap deflections, 
both positive and negative. 

It is to he observed that each critical curve is composed 
of three distinct parts: a substantially flat top and two 

steep sides. These three portions correspond to three differ- 
ent conditions on the airfoil determining the minimum pressure 
peak. 


The thin airfoils considered are characterized by having 
large additional velocities near the nose and but a moderate 
rise in the basic velocity distribution approaching the maxi- 
mum-velocity point. The combination results in a, velocity 
distribution which is double— peaked for lift coefficients 
differing more than a small amount from the design lift, 

The velocity peak at the nose appears suddenly, the after 
velocity peak still remaining; there is no continuous tran- 
sition. For a certain range of lift coefficient about the 
design lift the forward velocity peak is less then the rear 
•ne; for this range the critical Mach number is determined 
by the velocities near the maximum— velocity point of the base 
profile modified by the small additional velocities for this 
region. This latter velocity peak will be termed the mid- 
peak velocity to distinguish it from velocity peaks which 
appear at hinge points on flapped airfoils to be mentioned 
later. This ran^re of lift coefficient corresponds to the 
top of the critical curve. The small additional velocities 
produce a near linear change in Mach number with lift coef- 
ficient, giving to the top a slight slope as shown in figure 
2 . 


At some lift coefficient, the velocity at the forward 
pressure peak becomes equal to the veil city at the midpressure 
peak; then for large increments of lift from the design lift, 
the peak velocity at the nose is the maximum and changes 
rapidly with change in lift coefficient due to the large 
additional velocities. The result is that the critical Mach 
number falls rapidly with increasing increment in lift coef- 
ficient. The lift coefficient for which the velocities at 
the nose and midpressure peaks are eoual is the point of 
intersection of the top and the side; it is the point at 
which the absolute maximum velocity jumps from a position 
aft on the airfoil section to a point in the immediate 
vicinity of the nose and will be denoted as the "declination 
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point" on the critical curve. There are two such declination 
points, for positive and negative lift increments. For posi- 
tive lift increments the velocity peak at the nose appears on 
the upper surface; for negative lift increments it appears 
on the lower surface. 

Since airfoils with flaps are equivalent to the original 
airfoil with modified camber, the effect is to change the 
basic velocity distribution while leaving the additional 
velocity distribution the same. The critical curve, then, 
for the flapped airfoil is very similar to the unflapped 
but shifted as a whole by cj, The steep sides are 

D o 

practically parallel, for the incremental basic velocity at 
the nose is negligible. The only appreciable change occurs 
in the top and the declination points. 

For flap deflections small enough that the velocity peak 
at the hinge point is less than the midpeak velocity, the 
top of the critical curve will have essentially the same 
slope, since the airfoil has the same additional velocities, 
but is shifted up or down depending upon whether there is a 
negative or positive flap deflection. The increment of lift 
coefficient between the upper and lower declination points 
will consequently vary. 

For flap deflections large enough that the pressure peak 
at the hinge point is greater than the midpeak pressure, the 
top of the critical curve has a lesser slope than before, the 
additional velocities at the hinge point being less. It is 
to be noted that for negative flap deflections the peak 
velocity at the hinge point appears on the lower surface 
where the additional velocity increment with lift coefficient 
is negative; consequently, the top of the critical curve has 
a reversed slope to that of the unflapped and positively 
flapped airfoils. This is clearly shown in figure 2. 

The locus of the upper and lower critical Mach number 
declination points together with the top of the critical 
curve form what may be termed an "optimum critical curve," 

All points on this curve correspond to flap deflections for 
which, at a given lift coefficient, the airfoil section 
achieves the maximum possible critical Mach number, as can 
be seen in figure 2, 

The extension of the region of high critical Mach number 
by means of flaps is done in two parts. The first, where the 
hinge pressure peak is less than the midpeak pressure, is a 
near linear extension of the original curve. In this region 
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a reasonable gain in lift coefficient is realized without too 
much sacrifice in the critical Mach number. In the second 
part, after the break, the hinge peak predominates; the 
sacrifice in critical Mach number for increased lift is 
correspondingly greater. 

In figure 3 are shown the optimum critical curves for 
the NA.CA 64—, 65—, and 66— series of airfoils having thickness- 
chord ratios of 0,06, 0.08, 0.10, and 0.12 with 0.10c, 0.20c, 
and 0.30c plain flaps. The critical curves of the original 
airfoils without flaps (6=0) are in each instance shown as 
dotted lines. The rest of the critical curves for other flap 
deflections are not shown, as was done in figure 2, but rather 
only the locus of the declination points. These loci, as 
noted previously, form optimum critical curves for the given 
airfoils with flaps. Each of the curves in figure 3 has a 
form directly analogous to the optimum curve of figure 2, 
having the typical extension, for small flap deflections, 
of the h igh— cr it i cal— sp eed 1 if t— co ef f i ci ent range. 

There are two predominate variations to consider: section 

thickness and airfoil family. By family, reference is made 
to the position of the velocity peak on the base profile. 

The curves show a general trend to higher critical Mach numbers 
and smaller 1 if t— co ef f i c i ent ranges for the thinner sections. 
The extension of the 1 if t— coef f ic i ent range by means of flaps 
shows this same tendency though to a lesser degree. The 
decrease in the extension of the 1 if t— co ef f i c i en t range for 
the thinner sections is the result of higher velocity peaks 
at the hinge point, together with a smaller variation in 
velocity along the base profile so that the break in the 
curve occurs sooner. The larger peak velocities for the thin 
sections are due to the fact that the value of the peak 
velocity varies inversely as the section thickness at the 
hinge point. 

The additional velocity distributions of all the 6— series 
airfoils considered are substantially the same over the rear- 
ward 60 percent of the chord. For this reason, all the slopes 
•f the tops of the critical curves in any given airfoil family 
are essentially the same (cf. fig. 3(a) to 3(d)). The slope 
for different airfoil families will be different as indicated, 

( cf . fig. 3(d), 3(h), 3(l)) the slope decreasing for a rear- 
ward movement of the maximum velocity poixit on the base 
prof il e . 

For the different families a rearward movement of the 
maximum velocity peak is indicative of three changes in 
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section characteristics; a lower maximum velocity and flatter 
profile velocity, greater additional velocities at the nose, 
and a thicker section at the hinge point giving lower velocity 
peaks there. Thus, a rearward movement in pressure peak re- 
sults in slightly higher critical speeds with, however, a 
smaller range of lift coefficient in the high critical region 
together with a smaller extension of this range by the fla,ps. 

The extension of the critical curve of an airfoil by 
means of flaps will vary with the flap— chord ratio, the dif- 
ferent optimum critical curves being, however, quite similar. 

A study of the curves given in figure 3 indicates that, of 
the flaps considered, the 0.20c gives the best over— all re- 
sults . 

Experiments carried out in the Ames 1— by 3i— foot high- 
speed tunnel, on airfoils equipped with orifices for the 
determination of pressure distributions, have given some in- 
sight into the validity of theoretical calculations of 
critical Mach numbers. The sections considered had thickness- 
chord ratios equal to 0.12 and 0.15, the results being com- 
pared with those given in reference 1. It was foxmd that, 
throughout the high portion of the cr it i cal— sp eed curves for 
the low— drag airfoils, excellent agreement was obtained ex- 
cept that the points of declination of the curves were at 
lift coefficients somewhat beyond those predicted, especially 
for the case in which the peak pressure point liras forward on 
the upper surface. The source of this discrepancy is probably 
twofold: some error exists in the estimation of the true 

lift coefficient, and the Karman— Ts ien correction formula 
is erroneous in the immediate vicinity of the nose. It is 
possible to derive a modification of the Glauert correction 
for lift, by means of the Karman— Tsien pressure formula, and 
the change brought about in lift coefficient can be shown to 
be in the right direction to agree with experiment. This 
change, however , is small and somewhat laborious to apply and 
therefore was not used. That errors in predicted pressures 
should exist near the nose of an airfoil follows from the 
fact that the Karman— Tsien formula is postulated on the use 
of a pressure-density relation which holds for conditions not 
differing greatly from those in the free stream. Near the 
nose, in the vicinity of the stagnation point, some dis- 
crepancy would be expected to occur, and this is confirmed by 
the experimental pressure distributions. 

Considerations of theory and exper imental results both 
indicate that the theoretical results are conservative. It 
is to be expected, however, that, for airfoils with flaps. 
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the predicted extensions of the curves should he of the right 
order of magnitude. . _ 


CONCLUDING REMARKS 


Theoretical calculations show that the use of plain flaps 
on airfoil sections such as are considered in the present re- 
port will serve to increase appreciably the range of the high- 
critical Mach number region characteristic of low— drag air- 
foils. It is to be expected, judging from what experimental 
evidence is available, that the results will underestimate 
the extent of the critical curve lying between the declination 
points but the predicted extension, in this portion of the 
curve, should be of the right order of magnitude. 
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APPEND IX 
List of Symbols 

a Q velocity of sound in the free stream 

c chord length of airfoil 

Cj section lift coefficient 

° 0 I° w ~ s P ee d section 

c 7 incremental addit 

*a6 

c incremental basic 

M Mach number (velo 

Pj^ pressure coeffici 

PM=0 pressure coeffici 

Pf pressur e~coeff ici 

P chordwise lift di 

coefficients on upper and lower surface of airfoil) 

P a chordwise lift distribution produced by additional load 

distribution 

P-jj chordwise lift distribution produced by basic camber- 

line loading 

Pg incremental load distribution produced by flap deflection 

P a g incremental additional load 

flap delfection 


lift coefficient 
ional lift coefficient 
lift coefficient 

city divided by velocity of sound) 



ent under low— speed conditions 

ent d is tr ib\i.t ion over base profile 

stribution (difference between pressure 


distribution produced by 
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P-bg incremental basic load distribution produced by flap 
deflect ion 

basic load distribution for infinitesimally thin airfoil 

0 P ^ 6 incremental basic load distribution for infinitesimally 
thin airfoil 

p static pressure 

q dynamic pressure (^~pV 8 ^ 

r radius of curvature 

Au increment of local velocity produced by basic lift 
distribution on airfoil 

v local velocity 

Av a increment of local velocity produced by additional lift 
distribution on airfoil 

V 0 velocity of the free stream 

x distance along chord measured from leading edge of 

airfoil 

y ordinate of camber line measured from chord line 

Y ratio of specific heats (Cp/c v = 1.4) 

6 angle of flap deflection 

9 variable defined by equation 2x = c( 1 — cos 9) 

A tangent of angle 8 

p dens ity 

Subs cr ipt s 

o free— stream conditions 

cr critical conditions 

L lower surface of the airfoil 


U 


upper surface of the airfoil 
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Figure 1.- Equivalent camber line for airfoil with plain flap. 
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Figure 2.- Critical Mach number M C x -variation with high-speed lift 
.coefficient c Z j| for various flap deflections 6 NACA 

651-210 airfoil section, .20c plain flap. 
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Figs. 3a, b 
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Figure 3.- Optimum critical Mach number M« x variation with high-speed 
section lift coefficient c>y for .10o, .20c and .30c plain 
flaps, (a) NACA 64i-206 airfoil section. 
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(b) NACA 64 1 -208 airfoil section. 

Figure 3.- Continued. 
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(f) NACA 65i-308 airfoil section? 

Figure 3.- Continued. 
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(i) NACA 661-206 airfoil section. 

Figure 3.- Continued. 
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Figs. 3k ,1 
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airfoil Bection. 

Figure 3.- Continued. 
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(l) NACA 66i-212 airfoil section. 

Figure 3.- Concluded 













